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In contrast to the coupling parameter in the usual perturbative QCD (pQCD), the coupling 
parameter in the analytic QCD models has cuts only on the negative semiaxis of the Q 2 -plane 
(where q 2 = — Q 2 is the momentum squared), thus reflecting correctly the analytic structure of the 
1 spacelike observables. The Minimal Analytic model (MA, named also APT) of Shirkov and Solovtsov 

I , removes the nonphysical cut (at positive Q 2 ) of the usual pQCD coupling and keeps the pQCD cut 

■ discontinuity of the coupling at negative Q 2 unchanged. In order to evaluate in MA the physical 

" QCD quantities whose perturbation expansion involves noninteger powers of the pQCD coupling, 

a specific method of construction of MA analogs of noninteger pQCD powers was developed by 
Bakulev, Mikhailov and Stefanis (BMS). We present a construction, applicable now in any analytic 
QCD model, of analytic analogs of noninteger pQCD powers; this method generalizes the BMS 
approach obtained in the framework of MA. We need to know only the discontinuity function of the 
analytic coupling (the analog of the pQCD coupling) along its cut in order to obtain the analytic 
analogs of the noninteger powers of the pQCD coupling, as well as their timelike (Minkowskian) 
counterparts. As an illustration, we apply the method to the evaluation of the width for the Higgs 
decay into bb pair. 

43' 

O ; PACS numbers: 12.38.Cy, 12.38.Aw,12.40.Vv 

I. INTRODUCTION 

<rv 

It is well known that the perturbative approach to QCD (pQCD), while working well in evaluation of physical 
quantities at high momentum transfer (\q 2 \ > 10 1 GeV 2 ), becomes increasingly unreliable at low momenta (|g 2 | ~ 
\ 1 GeV 2 ). One of the main reasons for this is the singularity structure of the pQCD coupling parameter o pt (Q 2 ) = 
^(Q 2 )/ 71 " a t spacelike low momenta q: (0 <) Q 2 = —q 2 ~ 1 GeV 2 . This singularity structure does not reflect correctly 
the analyticity structure of the (to be evaluated) spacelike observables F(Q 2 ). The latter, by the general principles 
of the (local) quantum field theory [l|, Q , must be analytic functions in the entire Q 2 plane except on the cut on 
the negative semiaxis: Q 2 € C\(— oo,0]. Qualitatively the same analytic properties should have also the coupling 
parameter A\(Q 2 ) that is used (instead of a p t(Q 2 )) to evaluate the spacelike observables ^(Q 2 ). 

The first such analytic version was constructed in where the discontinuity function of pQCD (a) = 

Ima pt (Q 2 = —a — it) was kept unchanged on the entire negative axis in the Q 2 -plane. More specifically, the use of 
the Cauchy theorem for the (powers of the) pQCD coupling gives 
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X: 
S: 



„ 1 f°° da Ima"(-cr - ie) 

«pt(Q 2 ) = - / T^TyT, -, (1) 

where the integration is along the entire cut of the pQCD coupling in the Q 2 -plane (-co, +A 2 ), with < A 2 ~ 1 GeV 2 
being the (Landau) branching point, and 77 — ^ +0 (see fig. H]). Elimination of the unphysical (Landau) cut (0, +A 2 ) 



(integration path) 



2 2 
Q =0 A 2 \ Q -plane 



* The changes in comparison with arXiv:1106.4275v2: section II is extended (after Eq.(ll)); in section III: the part between eqs.(36) and 
(38) is new; appendix B is new; the part between eqs.(44) and (50) is new; the part between eq. (54) and the end of section III is new. 
^Electronic address: gorazd. cvetic@usm.cl 
tElcctronic address: kotikov@theor.jinr.ru 
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FIG. 1: The path of integration for the integral {TJ. 

in the above dispersion relation leads to the aforementioned Minimal Analytic (MA) 1 coupling 

„,,, „ 1 f°° da Ima^ti—a — ie) 
4 MA) W 2 )- t^t^ " (n = l,2,.... (2) 

It is named also Analytic Perturbation Theory (APT), It is applied usually in the MS renormalization scheme, 

or in truncated versions of that scheme. The method of equation (J2]) allows us to evaluate the MA coupling analogs 
A^ A) {Q 2 ) of the pQCD coupling powers a™ t (Q 2 ) even when n is noninteger (n4 v). 

MA gets its only free parameter, the QCD scale A, fixed by the requirement that it reproduce high energy QCD 
quantities ( | Q 2 | ~ 10 1 GeV 2 ) , and in this regime it gives good results, @ . It gives good results for the Bjorken polarized 
sum rule (a spacelike quantity) even at low Q 2 , [lj, [Hj], although at very low Q 2 ~ 0.1 GeV 2 it apparently requires 
a modification T3j. Also due to duality violations we should expect that, at low a ~ 1 GeV 2 , the discontinuity 
function pi(cr) = \va.A\{— a — it) in analytic QCD models will deviate significantly from the pQCD counterpart 
p^' t \a) = Ima p t(— a— ie). Another reason for the need of such a deviation is the apparent inability of MA to reproduce 
the correct value of the well- measured (timelike) low-energy QCD observable r T , the strangeless semihadronic decay 
ratio of the r lepton. Its present-day experimental value is r T (exp.) = 0.203 ± 0.004, [14j, |l5(. 2 In MA, the predicted 
values are in the range of 0.13-0.14, |J,[lg|, unless the values of the current masses of the light quarks (m u , m d , m s ) are 
abandoned and effective quark masses ~ 0.25-0.45 GeV are used instead [ijj ■ This numerical loss in the size of r T in 
MA appears to be connected with the elimination of the unphysical (Euclidean) part of the branch cut contribution of 
perturbative QCD, while keeping the discontinuity along the rest of the cut unchanged 18]. Furthermore, perturbative 
QCD models which are simultaneously also analytic (anpQCD), have also been investigated, [l9j], and they turn out 
to give too low r T value (r T < 0.16) unless their beta-function is modified in such a manner as to give convergence 
only in the first four terms of expansion, followed by explosive growth in the subsequent terms due to a rather singular 
choice of the renormalization scheme. 

Therefore, in general analytic QCD models, we must expect the following form of the dispersion relation: 

MQ 2 )= 1 - [ + °° ^ 

with p\(<j) = ImAi(-a-ie) deviating from p < f t \a) at low er, and (0 <) M 2 hr ~ M 2 being a threshold mass of the cut 
(— oo, —M 2 hr ) of Ai(q 2 ) in the Q 2 -plane. Having in such general analytic QCD models no direct relation of a pt (Q 2 ) with 
Pi(a), the general method of obtaining A n (Q 2 ), the analytic analog of the power a pt (Q 2 ) n , is not as straightforward 
as in equation 1^ in the case of MA. In [13,[2l|, the higher power analogs A n (Q 2 ) for integer n's and for any analytic 
QCD were constructed as linear combinations of logarithmic derivatives Ak(Q 2 ) oc d k ~ 1 A\{Q 2 ) / d{h\ Q 2 ) fc_1 (k > n), 3 
such that the evaluation of observables in analytic QCD leads to suppressed dependence of the evaluated truncated 
analytic series when the number of the terms in the series is increased. For the construction of A n (Q 2 ) (and its 
timelike conterpart 2l„(cr)), only the knowledge of pi(cr) (or equivalently, of A\(Q 2 )) is needed. The construction of 
higher power analogs A n , not as powers of Ai but rather as linear operations on Ai, has an attractive functional 
feature: it is compatible with linear integral transformations (such as Fourier or Laplace) [23| . 

It turns out that some observables, in particular mass-dependent ones, have pQCD expansion which involves 
noninteger powers a p t(Q 2 ) u ■ In order to evaluate such observables in any analytic QCD, we need to construct their 
analytic analogs A V (Q 2 ) (and their timelike conterparts 2l„(cr)). In the case of MA, a method of calculating such 
quantities was developed and applied in j24l - [27| (for a review, see [IH), their method being different from the direct 

evaluation ^ with n v [y noninteger). The analytic properties of their couplings A^ A ^ (Q 2 ) can be seen more 
clearly than in the formulas @, but numerically they are equivalent. In the present paper, we present the method 
of construction of _4„((2 2 )'s that is applicable in any analytic QCD model. Below we will demonstrate that, within 



1 Another, somewhat different, approach, performs minimal analytization of dlna p t(Q 2 )/rflnQ 2 function, [fj. An analytization using 
Borel transform of observables and the minimal analytization approach, was presented in Q. For reviews of various types of analytic 
QCD models, see [UH]]. 

2 r T represents the QCD part of the decay ratio R T , i.e., (r T ) p t = a p t + C(a 2 t ). 

3 The relations between ^^'s and An's allowing for recurrent construction of «4 n , for integer n, were given also in fil. |2'J . within the 
context of the MA model of 
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the MA model of Shirkov and Solovtsov, our approach gives the same result as the approach o f I24M26I ] in the leading 
(one-loop) order of perturbation theory. Above the leading order in MA, the results of [24l - |26j | are represented as 
certain expansions via the leading order results. Such types of expansions are absent in our approach. 

In sections [TTI and HU we derive the spacelike analogs A U (Q 2 ) of nonintegcr pQCD powers a p t(Q 2 )" ', as functions 
of pi(cr) of the analytic QCD model. In sec. [ill the construction leads us first to (noninteger counterparts) of the 
logarithmic derivatives, a p t^(Q 2 ), and their analytic analogs A U (Q 2 ). In sec. Mil we relate a p t^(Q 2 ) with the pQCD 
(noninteger) powers a p t(Q 2 ) lJ+m ] this relation is derived in appendix IA"1 This allows us to evaluate the mentioned 
spacelike observables J-(Q 2 ) in analytic QCD using either the analytic analogs A V {Q 2 ) of a pt ,u{Q 2 ), or the analytic 
analogs A V {Q 2 ) of a p t(Q 2 ) u . Furthermore, we construct the timelike counterparts 21,, (cr) and QL v (a) of the spacelike 
couplings A U (Q 2 ) and A V {Q 2 ). In sec. |IV] we apply, as an illustration, the presented method to evaluation of a 
timelike quantity, the width of the Higgs decay into bb pair, T(H bb). The corresponding spacelike quantity has 
a perturbation expansion which involves noninteger powers of a p t, due to the 6-quark mass anomalous dimension. 
We present the results for T(H — > bb) as a function of the squared Higgs mass s = Mjy, for various analytic QCD 
scenarios and in pQCD. In sec. [V] we present conclusions. 



II. LOGARITHMIC NONINTEGER DERIVATIVES OF EUCLIDEAN COUPLING IN ANY ANALYTIC 

QCD MODEL 

As mentioned in the Introduction, we will start with the analytic Euclidean coupling A\{Q 2 ) and its logarithmic 
derivatives in a general analytic QCD model. Such a model is determined (characterized) fully by the discontinuity 
function 

pi(cr) = Im A 1 {Q 2 = -a - is) , (4) 

defined for a > 0. Usually, the discontinuity cut is nonzero below a threshold value — a < — M 2 hr where M t hr ~ M„. 
The application of the Cauchy theorem to the function A%(Q 2 )/(Q 2 — Q 2 ) in the complex Q 2 -plane, along the closed 
path made of a very large circle and two segments just below and above the cut, gives us the well known dispersion 
relation for the analytic Euclidean coupling A\ 



1 P + t 
Ai{Q 2 ) = - / day! 

n Jo (v 



Pi{°) 



+ Q 2 ) ' 



(5) 



where q 2 = —Q 2 is non-Minkowskian, i.e., Q 2 can have any value in the complex plane except the cut (— oo, — M 2 h] .}. 
The logarithmic derivatives are defined as 

where /?o is the first coefficient of the (3 function: (3 — (1/4)(11 — 2n//3); Q 2 da pt /dQ 2 — ~(3 a 2 t + 0(a pt ). We note 
that for n = equation (j6]) gives A\ = A\ . We can write the logarithmic derivatives in the following form: 

~ o 1 f°° do 1 d n ( z \ 

A ^ iQ ) = -J a V^ (g) OT n + l) d (ln^ \lT~z) —~ ■ W 



=<?/Q 2 



It turns out that the integrand is the known polylogarithm function 



777 / \ / 7 \ 71 OO 

a / z \ id 



= ( z T) E (-l) m+1 ^ m = E (-l) m+1 m n z m = (-l)U. n (-z) , (8) 



d(lnz) n V 1 + z J \ dz, 
which brings equation ([7| in the following form: 



1 (-1) f°° da _ ,_, TJ f 



A n+1 (Q 2 ) = - ; ' / — Pl (a)U_ n (-a/Q 2 ) . (9) 
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This relation is valid for n — 0, 1, 2, .... Analytic continuation in n >— > v gives us 4 the logarithmic noninteger derivatives 

A v+ i{Q 2 ) = [ £*(<r)LU (-£) (-1 < u) . (10) 

We note that the integral converges for v > — 1. Namely, at high a (|z| 3> 1 where z = cr/Q 2 ) we have in the 
integrand of equation fj 10[) : pi(cr) w '(cr) ~ In -2 cr ~ ln~ 2 z and Li_ 1/ (— z) ~ ln _1/ z (for noninteger ^). Therefore, 
the integral converges at cr — > oo if f > — 1. The integral obviously converges at low cr, too. 

In principle, a continuation to arbitrary v 1 as performed by the transition from equation (j9]) to equation (|10|). could 
in principle miss some terms, such as terms proportional to sin (vk). However, such terms will be excluded because 
they are finite oscillatory when v — ¥ ±oo. 

It is interesting that the recursive relation 

X +2 (Q 2 ) = , i'^w/^ A+ilQ 2 ) . (ii) 

p Q {v + 1) d\nQ z 

which for positive integer v = n = 0,1,2,... is a direct consequence of the definition ([TJ, remains valid even for 
noninteger v as a consequence of the relation (I10p and the known 5 relation z(d/ dz)lA- v (z) = Li_„_i(z). 

We can recast the result (|10[) into an alternative form involving the spacelike coupling A\ instead of the discontinuity 
function pi(cr). This can be performed in the following way. 

We can use the following integral form of Li_„ function ([I(|) 6 appearing in equation (fT0|) : 

The last expression on the right-hand side was obtained by the change of variable t = ln(l/£). Since we have in our 
result (fT0|) Li-ix with —\<v (and not just: — 1 < v < 0), we extend the integral representation to higher v > 0. 
This is achieved by using in equation (1121) the aforementioned relation (d/dln zVLi_„ = Li_„_i. We thus obtain, for 
v = n + 6, with < 5 < 1 and n = —1, 0, 1,2,.. ., the following integral form, |3l| : 



Li_„_,5(z) = 



d\n ; 



r(i-5)7 i-zc n U 



(n= -1,0,1,...; 0<5< 1) . (13) 



Inserting the representation (|13p . for f = rt + (5, into our general formula (|10[) . and exchanging the order of integration, 
gives us 

The last integral over da is the spacelike coupling ^li(Q 2 /^) due to the dispersion relation (J5)). Therefore, we obtain 
the alternative form of the result (fTU|) . for ^ = n + 5, with < 5 < 1 and n — —1, 0, 1,2,.. ., 



(16) 



^r(i + i/)r(i - 5) V dinQV Jo £ v / ; U 

1 r(l + (5) sin(^) x J '°° 



/3gT(n + l + 5) (ttcT) V dm <2 



where the last form (|16[) was obtained from the previous one by the substitution t — ln(l/£) and using the identity 
T{l + S)T{l-S) = -K5/sai{irS). 



4 In Mathematica |28)l , the Li_„(z) function is implemented as PolyLog[— u, z] . However, at large \z\ > 10 7 , PolyLogf— v, z] appears to 
be unstable. For such z we should use the identities relating Li_„(z) with Li_„(l/z), which can be found, for example, in |29l |. 

5 This relation can be obtained, for example, by applying d/dlnz to the power series Li„/(z) 

scr. — / , Z Tfl. 

6 Equation ||12|I can be proven by expanding the integrand in powers of e~ 4 and using the basic integral expression for the Y(u ) function 

(where v = —v > 0): du e~ u u v ~ 1 = T(v ). In this way, the (convergent for \z\ < 1) series 5Zm=i zm " m ~ 1 ' ls generated, which is 
just the polylogarithm function Li^/ (z) = Li-v(z). 
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Furthermore, we will now prove that the obtained result (fT0|) (equivalent to equations (IT5|) and (jTHJ) ) reproduces, 
in the specific case of the Minimal Analytic model (MA) of (JOf at one-loop level, the explicit result obtained in [24j 7 



Li_„(A7Q 2 



/3o +1 Vln" +1 (g a /S ) ty + l) 



(17) 



where the scale A appears in the one-loop MA analytic coupling Ai(Q 2 )( MA ' 1 ^ and in its discontinuity function 



2\(MA,l-£) _ J_ 



1 



A 2 



A) Vln(g 2 /A ) (Q 2 - A ) r 



Pl {at- l) = Tma pt (-a-ie)^- £ >=lmA 1 (-a-ie)^ 1 - e > 



-Im 



Po (ln((i/A) - itt) A) (ln 2 (a/A) + tt 2 ) 



(18) 



(19) 



When replacing _4i(<3 2 /£) in the integrand of the expression (fT5"j) by the second term of the expression (JTHJ) for 
Ai(Q 2 1 '^)( MA ' 1_£ ) J and using the integral form (fT3|) for Li_„, we obtain immediately 



-In 

o s 



1\ 1 



(-1)A Z 



i) /?o( Q 2 /c _A 2 ) /?o +1 r> + i) 



( " Li_ tf (A a /Q 3 ) . (20) 



On the other hand, when replacing A\ {Q 2 /(,) m the integrand of the expression (|16l) by the first term of the expression 
(dHJ) for Ai(Q 2 /£,)( MA ' 1 ~ l \ we obtain in a direct manner 8 



1 r(l + 6) sin(7r<5) 
%T{n+l + 6) (tt<5) 



rflnQ 2 



n+l 



(ft 1 



1 



t 4 A> 



i + ln(Q 2 /A 



^ +1 ln" +1 (Q 2 /A~) 



(21) 



Combining the results ([20| and (|2T|) . we obtain the full result ([T71) for A V +\(Q 2 ) in the one-loop approach of MA, for 
any noninteger v such that —\<v (when v is nonnegative integer, the limit 6 — > can be made in the derivation). 
This (one-loop MA) result, obtained for the first time by Bakulev, Mikhailov and Stefanis (BMS) in (24j, has several 
interesting properties, as pointed out in [25] (their equations (3.14)-(3.19)). The result (|T7)) is explicit and allows us 
to apply it even for v < — 1, and even for complex i>; this is a kind of analytic continuation in v. We can thus use 
this result, by adding and subtracting it from of our general integral expression (|10j) . thus extending the (/-regime of 
applicability of our expression 



X + i(Q 2 ) = A V+1 (Q 2 )^ A ^ + -- ( 1} 



tt /3£I> + 1) J 



da 



{!-€) 



LL 



-— ] (-2- ,,.) . (22. 



where .A„ + i(Q 2 )( MA ' 1_£ ) and p\(o) v ^ ' are given in equations ([TT]) and respectively. Now the integral converges 
also for — 2 < v < —1, because, due to asymptotic freedom, the difference \fii(cr) — Pi{&)pt~ ] behaves at large a as 
~ lnlncr/ln 3 (T and not as l/ln 2 cr. Further, the expression (122(1 implies that Ao(Q 2 ) [= lirrv_>_i.4. I /+i {Q 2 )] = 1 for 
all complex Q 2 , because: U- V (z)/T(v + 1) -> when v -> -1, and A(Q 2 ) (MA ' 1- ^ = 1- 



7 Note that a in [244271 ] corresponds to our floa (with our fit) = (1/4)(11 — 2n.j-/3)); their A v +i corresponds to our /3q +1 A v +i ; they use 
the transcendental Lerch function notation z <3?(z,i/ ,1) = F(z,v ) for the polylogarithm function Li / (z). On the other hand, An in 

corresponds to analytic analogs of a™ = 7r"a n , i.e., their A n corresponds to our Tr n A n - 

8 We can use the integration variable y = t/to, where to = ln(Q 2 /A ), and the exact solution of the following integral: 

dy it 



L 



y S (y + l) sin(7r<5) 



where < 5 < 1 
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III. ANALYTIZATION PROCEDURE FOR OBSERVABLES WITH NONINTEGER POWERS OF 

COUPLING 

In QCD we encounter often spacelike (Euclidean) observables T(Q 2 ) whose perturbative expansion starts with a 
noninteger power a p ° t 

^(Q 2 ) P t = a pt (Q 2 r + F iapt (Q 2 r +1 + J- 2 a pt (Q 2 r +2 + • • • (23) 
The general analytization procedure of the pQCD-evaluatcd observables with integer powers is 

apt,n+i h> A n +i (n = 0, 1, 2, . . .) , (24) 
where a p t, n +i are the logarithmic derivatives of the pQCD coupling a p t 



and 9 /3o is the first coefficient of the /3-function 

dapt{ ^ 2) - a(a pt ) = -A,«4t-^-A4t-ft<4t-^oS 



= -A)»pt (1 + CiOp t + c 2 a 2 t + cga^ + c 4 a* t + . . . ) fc,- = ^ 



(26) 



In the case of observables whose pQCD-evaluated expressions are the (truncated) expansions equation ([23]) with 
noninteger Uo, the analytization procedure (IM|) is naturally extended to (|20l. [2lj) 

2pt,w+i ^ j ( 27 ) 

where the expression for A v +\ is given in equation (|10p . Therefore, at this stage, the problem of evaluation of such 
observables in anQCD is reduced to re-expressing the noninteger powers a pt in pQCD expansion (|23|) in terms of the 
logarithmic noninteger derivatives 'a P t,v+m{Q 2 ), in order to perform the subsequent analytization via equation (|27[). 
Stated otherwise, we find first the coefficients k m (v) of the relations 

oo 

«pt,, = o£t + E M")C" - (28) 

m— 1 

and, as a consequence, the coefficients k m (v) of the inverse relations 

oo 

(29) 

m— 1 

The expressions for the coefficients k m {v) (and k m (v)) are derived in appendix lAl see equations (|A6|) - (|A9|) . (|A10[) and 
(IA11[) there for explicit expressions. There, the coefficients k m (n) and k m (n), for n integer, are obtained by solving 
the difference (recursion) equations relating k m (n + l) with k m (n), k m -i(n), etc. The solution for k m (n) (and k m {n)) 
is obtained in a form involving combinations of Gamma functions T(x) and their derivatives (up to m derivatives), 
at the values of the argument x = 1 and x = +n + ml (for m! = 1, . . . , to). In the obtained expressions, the integer n 
is then replaced by an arbitrary noninteger v (n i— > v). The latter step is an analytic continuation similar to the step 
ri H> v from equation (|9]) to equation (H~0|) . 



9 Naively, one might suppose that the analytization procedure, in the evaluation of observables .F(Q 2 ) = T>(Q 2 ) with integer powers of 
a p t, in any given anQCD model would be a™^ 1 >— > It turns out that, in those anQCD models whose Ai(Q 2 ) at high Q 2 differs 

from a p t(Q 2 ) by negative powers of Q 2 (~ (A 2 /Q 2 ) fc ), such naive analytization procedure leads to strong renormalization scheme (RS) 
dependence of the truncated (modified) analytic series T>( N ' (Q 2 )( m ) an , due to the contributions of power terms ~ (A 2 /Q 2 ) m to the 
derivative dT>W (Q 2 ) (m)an /aRS, see [p. 
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The relations (|29|) allow us to reexpress the expansion (|23|) in terms of a p t,i/'s 



'nipt 



(Q 2 ) + ^Fia p t^ +i{Q 2 ) + -T^ipt ,1/0+2 (Q 2 



(30) 



where 'mpt' stands for "modified perturbation series" and the coefficients T n are related with the coefficients T n of 
the original perturbation series (f2"3"f via relations involving the coefficients k m (i/q + n) appearing in the relations (j2T))) 



T\ = T\ + ki(vo), fi = Ti + T\ki(yv + 1) + k 2 {vo) , etc. 



(31) 



See equations (|A22[) - (| A25|) in appendix lAl for more relations involving higher orders. 

At this stage we apply the analytization procedure (f2~7|) to obtain the "modified analytic" (man) series 10 for the 
(dimensionless) spacelike quantity F(Q 2 ) 



•F(Q 2 )ma„ = A Uo {Q 2 ) + JiA 0+ i(Q 2 ) + -F 2 A +2(Q 2 ) 
l2\ 7 /n2^ 



(32) 



li 



where the expressions for A Vo +\{Q ), A Vo +2{Q ) are given, in any given anQCD, by equation (fT0|) 

On the other hand, an observable T(cr) that is related with a spacelike observable J r (Q 2 ) via the integral transfor- 
mation 



F{Q 2 ) = Q 2 

is timelike (Minkowskian). The inverse transformation is 

T "> - 55 



da T(cr) 

n (a + Q 2 ) 2 



'2 



Q 



'ttHQ ) > 



(33) 



(34) 



where the integration contour is in the complex Q 2 -plane encircling the singularities of the integrand, e.g., path C\ 
or C2 of fig. [2] Application of the tranformation (l34|) to the (modified) analytic series (|32|) then gives for the timclike 









. Q -pi 


/' 

/' 


\ 




i 


l-o-ie 


' ! — ' " 







FIG. 2: Paths Ci and C2 in the complex Q 2 -planc. 

quantity T(cr) the following (Minkowskian) "modified analytic series" 

7» man = 8^(0-) +^1^+1 (<J) +^2^+2^) + 



(35) 



The word "modified" is used here because we are analytizing the logarithmic (noninteger) derivatives of a p t [equation (125 M and not the 
(noninteger) powers of a pt . Our method leading to the expression d 1 1) makes the described former "modified" analytization approach 
more direct than the latter (equivalent) analytization approach involving (noninteger) powers [cf. equation 1 1521) later in this section.] 

Similarly as was argued in |2 j| in the case when vq is integer (vq = 1), it can be shown that the truncated series T{Q 2 ; fJ?)ma U , whose 
last included term is ~ ■A^+jvCm 2 ) an d the renormalization scale fj, 2 is used, has a systematically suppressed renormalization scale 
dependence when the order index N increases dJ r (Q 2 ; fJ. 2 )\^J n /d In fi 2 = 0(A vg +N+i)- This is really a systematic suppression, because 
in analytic QCD models we have the hierarchy |.4„ (/i 2 )| > |.4„ + i(/i 2 )| > ■ ■ • , for all complex fi 2 outside the cut. 
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where the timelike (Minkowskian) couplings 2l„+i(cr) are defined as 

^ +1 (a) = — -^ Au+1 (Q 2 ), (36) 

^1 J — a —ie H 

and the inverse transformation is 

x^sn-* £ ■ (37) 

Here we would like to note that the reexpression of the expansion ([25)1 to the one of equation ([50]) is a well-defined 
operation for quite convergent series. It is well known that mostly the QCD series are assumed to be asymptotic (see, 
for example, [32j]). and some arguments to justify such a reexpression should be done, 12 despite the fact that we use in 
our analysis only the first several terms in the expansions (I25"|) and . Some cases of different types of asymptotics 
were recently considered in |33T ]. 

To show the correctness of the reexpression and, at the same time, to avoid an additional increase of the volume of 
the main text part of our paper, we consider in appendix [B] the standard Lipatov-type behavior [34j for the nth term 
of the expansion (|23|) and recover the similar behavior for the nth term of the expansion (|30|) . We show that there is 
even a slight weakening of the rise for the nth term of the expansion Q30p in comparison with the original one in the 
expansion (l23l) if the coefficients in (|23|) are nonalternating in sign. 13 

The parameter a p t^ + n {Q 2 ) of the expansion (|3"0)) is very close to the corresponding one a l £ t +n {Q 2 ) in the high- 
momentum regime (i.e., when a p t-values are small). This observation is consistent with the similarity of the coefficients 
T n and T n shown in appendix [B] and another argument for acceptability of the reexpression of the expansion (|23[) to 
the one of equation (|30]) . 

We note also that, after the analytization procedure (|27[) . there is in general a significant supppression of the 
new parameters A Vo + n (Q 2 ) in comparison with a p t l! y +n(Q 2 ) and a,p° t +n (Q 2 ) (see a recent review and discussions 
therein) . 



A. Timelike (Minkowskian) coupling parameter 

Direct use of the expression (flO|) in the integral (|36|) gives us a double integral. When — 1 < v < 1, the order of 
integration can be exchanged because the resulting double integral is convergent, and we obtain 

f-ll r°° da' r-a+ie ,q'2 

**>+M = 9 -irv \n / ^tp^') / ^s-LM-oVQ 2 ) • (38) 

2m(3»T(v + 1) J a' J-a-ie Q 

Now we can use the fact that \A- V (z) is a function with a branch cut discontinuity [1, +oo) in the complex z-plane, 
and the earlier mentioned relation Li_„(z) = cHAi— v (z) / 'din z , to obtain the identity 

^-l — Li-„(z) = -ImLii_„(«-ie) = - — ^ . (In k)~ v Q[k - 1) , (39) 

where on the right-hand side is the Heaviside step function, and the integration on the left-hand side is along the 
contour of radius \z\ — k over the angles $ = arg(z) from +0 to (2tt — 0), see fig. [3] Using the relation (|39[) in 
the integration over Q 2 on the right-hand side of equation (|38|l . we obtain the simplified expression for the general 
Minkowskian coupling 2l„ + i(cr) in any analytic QCD and for any real v in the interval — 1 < v < 1 

~ . „ sin(7w) f 00 dw . . . 

^+iW = ^r — Pi^e"), (-Ki/<1), (40) 

and we used here a new variable w = \n(a' /a). The integral (|4U|) is clearly convergent at w — > 0. It is also convergent 
at w — > +oo, because there p\(ae w ) ps / o^ pt ^(cre ! ") ~ 1/w 2 . 



We thank to anonimous Referee who drew our attention to this possible problem. 

A decrease of T n in comparison with T n for n = 3 and 4 has been earlier observed also in [3f| . 
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z-plane 








> K+ie 




/ 

/ 







FIG. 3: Integration contour for the integral (1391) . 



The case of 21 1 (a) is obtained in the limit v — » +0 of the above expression 



2li(a) =2ti(a) 



1 



<iu> pi(<re u ') 



1 



do-' , 



(41) 



which is a well known result. For Minkowskian couplings 2l„ + i = ^s+i+nW) with higher index f — n + 5 (n 
0, 1, 2, . . .; and < 8 < 1), we can use the recursion formulas 



(-1) 



-21„+i(ct) , 



/3 (^ + 1) din. a 

which can be obtained from the relations (jTTJ) and (|36l) , and obtain (see appendix [C] for derivation) 



sin(7r(J + n)) f 00 dw 
Tr 2 (8 + n)p s +n J 



pi(ae w ) - Pl {a) 



w dp\{a) 
il diner 



,,n— 1 jtl—\ 



(n- 1)! d(ln cr)"- 1 



(42) 



(43) 



where ^ = n + 5, with < <5 < 1 and n = 0, 1, 2, When n = 0, the expression in the brackets in equation (|4"3"f is 

[ / oi(cre'"')], and <5 = ^ in this case varies in a larger interval — 1 < S < 1, i.e., equation (|4T))) . The integral (|4U1) is clearly 
convergent at u; — > +oo. It is also convergent at w — > 0, because the expression in backets behaves as ~ w n there. 

The version of this formula when 5 = is obtained by repeated application of the recursion formula (|4"2")l to the 
expression (l4"TT) 



o} n = ("!)" = 1 d"-Vi(a-) 

#?n! <9(lnfj)« fljn! Trd(lna)"™ 1 ' 



(n=l,2,...) 



(44) 



As we did in the previous section for the spacelike coupling A v +\, we derive now from our general timelike coupling 
43]) the explicit result obtained in [25[ for the one-loop MA case. First we rewrite the integrand in equation (|43|) 



n-l 

piK)-E 

fc=0 



w d pi(a) 



E 



w k d k p\{a) 



k\ d(ln cr) fe ^ A;! d(ln cr) fc ' 
Using the one-loop MA expression for p\{a) ( ^ t *'\ equation (p~9|) . we can represent (j45|) in this case as 

^dVi^)^ i T (-wr ^ {-wf 

— Im- 



(45) 



E 



tE- 



fe! d(lna) fe (ln(o/A 2 ) - i7r)"+! ^ (ln(a/A^) - ivr)^ 

1 



— Im- 



{-w) n 



(46) 



#o (ln(cr/A 2 ) - itt)" (ln(cr/A 2 ) + w - in) 
Putting the result (|46l) into the integrand in equation (1431) and performing the integration over w (again using the 
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integral given in footnote (J5J), we obtain 



1 



(5 + n)Ti(3 s + n+l 
1 



Im- 



ln(cr/A ) — in 



S+n 



(47) 



1 



ln 2 (<r/A') +7T 2 
The expression 



I//2 



sin 



artan(7r/ln(cr/A 2 ))^) . (48) 



where v = n. + S, with < <5 < 1 and n = 0, 1, 2, . . .. The expression (|48p is explicit and is, by (analytic in v) 
continuation valid for any v. It coincides with the result of BMS in 25] 14 and has, therefore, several interesting 
properties derived and specified on 25] (their equations (3.12)-(3.13) and (3.16)-(3.19)). Similarly as we did at the 
end of sec. IIIII for the spacelike coupling, we can now use this explicit MA one-loop timelikc coupling expression (valid 
now for any v) in order to extend the i/-regime of applicability of the general anQCD time-like coupling formula (I40j) 
[or: (|4"3"1) with n = 0] down to v w — 2 



(MA,l-£) 



JO 



dw; r 
it; 



Pl {ae w ) - Pl {<ye w ) 



{i-t) 
P t 



-2 < v < 1) 



where 2l„ + i(o-) 



(MA,l-() ■ 



is given in equation (|48j) . and /^(c)^. ^ in equation (fT9|) . We can apply the limit 



(49) 



-1 in 



equation (gSJ), and obtain 2t (cr) = 1 (for all a > 0), because Slo^)^ 1 "^ = 1. 



B. General form for the spacelike and the timelike observables 



We wish to stress that the formulas (fTU)) , (f2"2")l and (|4H|) , (|4*3")l , (l4*9"j) allow us to calculate the corresponding couplings 
A„+i and 2lj,+i for any real v > —2 and in any analytic QCD theory in which we know the discontinuity function 
pi(cr) = lmAi(Q 2 = —a — ie) (or equivalently: the coupling function A\(Q 2 )). 

We can define the combinations of Av+ n (Q 2 ys 15 which are analogous to the pQCD relations (|29|) under the corre- 
spondence (|2T|) 

A v = A u +^2 k m {v)A u+m (i/>-2). (50) 

TO>1 

As expected, it is easy to check that the analytic series (|32[) can then be rewritten in the form 

HQ 2 )™n = HQ 2 )an = A„ (Q 2 ) + T 1 A„ 0+1 {Q 2 ) + T 2 A„ 0+2 (Q 2 ) + --- ■ (51) 

Therefore, the comparison with the original perturbation series in powers of a p t(Q 2 ), equation (|23[) . gives us the 
correspondence between the pQCD and anQCD quantities 

a^t 1 -> A v+ i (52) 
for any real, in general noninteger, v > —2. Using the same combinations for the timelike couplings 21^+n's 

%, = %,+ ^k m {v)% v+m (v>-2), (53) 

m>l 

we can rewrite the associated timelike observable T(cr) of equation (|34jl in a form similar to the expansion (|35j) but 
involving the original J-j coefficients instead of JFj 

T{a) malL = T{a) &Il = % V0 {(T)+F 1 % Vo+1 {a)+F2^u o +2{(T) + --- . (54) 

In equations (|5"Tj) and (|54|) . the subscript 'an' now stands for "analytic". 



Their 2t„+i(o-) is our f3"+ 1( A u+1 (a), where ft = (1/4)(11 - 2n//3). 

We recall that the latter are defined via the integrals of equation d 1 1) or 1122 I I. 
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Equation (|50|) further implies that Aq(Q 2 ) — 1 (for all complex Q 2 ), because: (a) Ao(Q 2 ) = 1 as shown at the end 
of sectionHJ and (b) k m (0) = 1 as seen from the results (|A17j) - (|A20|) . Analogously, equation (|53|) implies 2lo(cr) = 1 
(for a > 0), since 2lo(< 7 ) = 1 as shown at the end of subsection IIII Al This means that Aq = 1 and 2lo = 1 in any 
anQCD, and this represents a consistency check of our method of construction of A u and % u . 16 

The results of our method allow us to obtain also analytization of powers combined with logarithms of the coupling 



a fc «pt(Q 2 ) 

'ptW J 111 "ptW J — Q^k 



a» pt (Q*)ln K a pt (Q*) - 

A v ,k{Q 2 ) \= (a v {Q 2 )\n k a{Q 2 )\ ]=^^2 (y > -2; k = 0,1,2,...) (55) 

L V / anJ OV K 

The right-hand side of equation ([55]) follows from the analytization rule (|52|) . applied separately to each power in the 
expression (a p ^ s — a pt )/8 (where 8 — > 0) when k = 1; when k > 2, the principle is the same. The derivative d k jdv k 
on the right-hand side of (|55|) is applied to each term on the right-hand side of the sum (|50j) . where we have to take 
into account that the ^-dependence is in the coefficients k m (y) and in the couplings A v + m (.Q 2 ) whose expression is 
given in equation (|22|) . 



IV. APPLICATION TO THE HIGGS DECAY WIDTH 



Following references [2J-|26(, in this section we apply the presented approach to the evaluation of the decay width 
of the (Standard Model) Higgs into heavy quark-antiquark (bb) pair: T(H — > bb) 

AT f~l 

T(H^bb)(s) = ^V~sT(s), (56) 

where Gf is the Fermi coupling constant, s = Mfj is the square of the Higgs mass, and T(s) is the imaginary part 
Imil(— s — ze)/(67rs) of the correlator of the scalar current Jb = m^bb 

U(Q 2 ) = i(4tt) 2 J dxexp(iqx)(0\T[J b (x)J b {0)}\0) , (57) 

where Q 2 = —q 2 , cf. (36l[37j. Later on in this section, we will see that the perturbation expansion of the corresponding 
spacelike quantity involves noninteger powers of a p t, due to the b quark mass anomalous dimension. Using the notations 
of [Hj], we can write the timelike quantity T(s) as a perturbation expansion 

T(s)=ml(s) fi + E*« a Pt( fi )j ' ( 58 ) 

where the square of the (spacelike) renormalization scale /i was chosen to be fj, 2 = s, and r?Tb(/i 2 ) is the MS running 
mass of the b quark. The corresponding spacelike quantity F(Q 2 ) is 



F(Q 2 ) = Q 2 f°° ^ 
Jo (°- 



da T(cr) 



Q 



2\2 



(59) 



and its expansion is written as 



F(Q 2 ) = m 2 {Q 2 ) M + J2 /™«pt(Q 2 )J , (60) 

Relations between the (dimensionless) coefficients fj and tj are given in 38] . 

The idea is to evaluate first the spacelike quantity F(Q 2 ), and obtain the timelike quantity T(s) (and thus the 
decay width) by application of the integral tranformation inverse to (l59l) [cf. also equations (|33|) - (j34)) ] 



We thank S.V. Mikhailov for pointing this out. 
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A. Running mass 



For this, we will use, in the expression (|60p . for the square of the running mass an expansion in (noninteger) powers 
pt( 



of a pt (Q 2 ). We recall that the renormalization group equation (RGE) for the squared MS running mass is 



m 7m(apt) = -ma pt ( 1 +5^7j«pt | > ( 62 ) 



din /i 



where the coefficients 7- (j = 1,2,3) of the mass anomalous dimension are known ([39|-|4l|]); for nt = 5, which applies 
in the case of the considered decay, we have: 71 = 3.51389, 72 = 7.41986, 73 = 11.0343. The 5-loop coefficient 74 
has not yet been calculated. Nonetheless, application of Pade approximants to the quark mass anomalous dimension 
7 m (flp t ) for rif = 5 indicates that 74 ps 12., and we will use this value. 17 Furthermore, the ft (j = 0, 1, 2, 3) coefficients 
(in the MS scheme) of the RGE (|26|) for the coupling a p t(Q 2 ) have been calculated explicitly, |42Tl45j |. and for nj = 5 
their values are: ft = 1.91667, ft = 2.41667, ft = 2.8 2 6 68, ft = 18.85 22. The 5-loop beta coefficient has been 
estimated in [46[ by Pade-related methods, and for rif — 5 the estimated value is ft = 165.161, which we will use 
here. 

Integration of the RGE's (|26p and (l6"2"j) gives for the squared running mass the solution 

m\{p?) = ml a;i(^ 2 ) (l + £ ^4tO* 2 )j (63) 

where m 2 is a renormalization scale invariant mass, v$ = 2/ft = 1.04348, and the coefficients Mj (j — 1,2,3,4) are 
functions of ft, Ck = ft/ft and 7^ (fc < j), and they are given in appendix [Pi For n/ = 5 these coefficients are: 
Mi = 2.35098; M 2 = 4.38319; M 3 = 3.87308; M 4 = -22.2155. The mass quantity m b is RG-invariant and can be 
obtained with high precision in the following way. The world average value of the QCD coupling parameter (in MS 
scheme) is a p t(M§) = 0.1184/7T [47j. The value of the running mass at its own renormalization scale, to^to 2 ) can be 
extracted from the heavy quarkonium physics. We will take the (central) value obtained in 48]: mb(m 2 ; nj = 4) = 4.24 
GeV. If taking the threshold between nf = 4 and n./ = 5 at fj, = 4.24 GeV, the aforementioned mass value and the 
world average value a{M%) — 0.1184/7T [43] lead to the values 18 at n/ = 5 

m b (ml;n f = 5) = 4.232 GeV , a pt (m 2 ; n f = 5) = 0.22542/vr . (64) 

Using these values in the relation (|63|) . we obtain the scale invariant mass 

m b = 15.330 GeV . (65) 



On the other hand, the values of the coefficients fj of the expansion (f60|) for j = 1, 2, 3 were obtained in [50|, and for 
j = 4 in [5l| (denoted as ft there): fi = 5.66667; f 2 = 51.5668 - 1.90696n /; f 3 = 648.709 - 63.7418n/ + 0.929133n 2 ; 
U = 9470.76- 1454.2871/ + 54.7826nj - 0.453744n^. The values for the here relevant case n f = 5 are: fi = 5.66667; 
h = 42.032; f 3 = 353.229; / 4 = 3512.2; respectively. 



B. Higgs decay 

We can now define the dimensionless ("reduced") spacelike quantity by dividing the expression F(Q 2 ) [equations 
59}, (|60j) ] by the RG-invariant scale m 2 , and using the expansion (|63"1) 



F{Q 2 ) = = a(Q 2 r + V T n a»°+ n (Q 2 ) , (66) 



17 Namely, applying to 7 m (a p t) (at nf = 5) the Pade approximants [3/l](a), [2/2](a), [l/3](a), and reexpanding in powers of a p t up to 
a^ t , gives us 74 = 16.4,9.2,10.2, respectively; the arithmetic average is 11.9, i.e., approximately 12. If repeating the same procedure at 
one order lower, we obtain from Pade approximants [2/1] and [1/2] the values 73 = 15.7,8.8, respectively, the average being 12.2 which 
compares favorably with the exact value 73 = 11.0343. The latter test give s us reason to except that 74 = 12. is a reasonable estimate. 

18 The discontinuity in the mass value at threshold can be obtained from 14911 . 
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where the coefficients F n are now the corresponding combinations of the coefficients fj and Mk 

?n = U + fn-lMx + ■■■ fxM n -X + M n • (67) 

For n f = 5 this gives: Fx = 8.01764; F 2 = 59.7374; F 3 = 480.756; Fa, = 4526.6. The expression F(Q 2 ) of equation 
(|66p is now the expansion of a spacelike quantity in noninteger powers (with vq = 2//?o = 1.04348) considered in the 
previous section, cf. equations (|2"3"j) . (|3U1) . (|3"2"j). The corresponding timelike quantity T(s) is 

Us) , ^ = ,/ (g ^ )(a) ^ (68) 



i2iV c G F V^/(4^V2) 



where s = Mj^ and we used the relation (|56[) . This quantity is then evaluated by the formula (|35|) . with ,F n 's 
(n = 1,2,3,4) determined by the relations pip , as explained in the previous section; or, equivalently, evaluated by 
the formula (|54[) . The evaluation can be performed in any analytic QCD theory, and even in perturbative QCD, 
simply by using in expressions (|40p and (|43[) the discontinuity function of the theory pi(cr) = InxAi(— cr — ie) (in 
anQCD) or p\{<j) — Ima p t(— cr — ie) (in pQCD). 



C. Numerical calculations 



For numerical illustration of our approach, we will consider here the discontinuity function pi(cr) to originate: (a) 
from the Minimal Analytic (MA) model of Shirkov and Solovtsov (also known as Analytic Perturbation Theory 
- APT); (b) the models which have, at high a > M^, the same pi(cr) as the perturbative QCD - this includes analytic 
QCD models of the type [52], 53] , and the perturbative QCD itself. We could contruct, in principle, such discontinuity 
functions by numerically integrating the RGE for a p t(Q 2 ) (in MS renormalization scheme and with an initial condition 
at Q 2 = Mg) over the complex plane of Q 2 and evaluating the imaginary part over the negative semiaxis. However, 
such an approach is cumbersome. We calculate Pi(cr) by evaluating a v t(Q 2 ) for complex Q 2 as a sum of the exact 
two-loop solutions a pt (Q 2 ,2 — £.) (which involve Lambert function, cf. |54l. |55|) as described in (5|| 



i P t(Q 2 ) = a P t(Q 2 , 2 - L) + Y,Cj^ t (Q 2 , 2 - I.) 



3=3 



where 



C 3 = c 2 , C 4 = -c 3 , C 5 = 



— C Cl C-i H Ca 

3 2 6 3 



i {-cxc% + c\c z - 2cxc 4 ) + 2c 2 c 3 + ic 5 



(69) 



(70) 



We truncate the series at j = 6. The last coefficient Cq depends also on C5 = P5/P0, which we do not know in MS 
scheme (even the estimates are not reliable), so we set C5 = 0. Since the Lambert function can be called upon in 
various numerical softwares, including Mathematica [28| . this high precision evaluation of p\(cr) = Ima pt (— a — ie) is 
fast. The two-loop coupling in terms of the Lambert function W±x, @, [55|, is 



(71) 



where Q 2 = \Q 2 \ exp(ic, 



ci [l + W Tl (z)} ' 

the upper subscript refers to the case < <j> < +tt, the lower subscript to — tt < <f> < 0, and 



1 

cie 





-Po/ci 


■^1 




exp 


—1 — q> 



(72) 



p t(M§) = 0.1184/7T [the corresponding usual MS scale (at nj — 5) is A = 0.213 GeV]. These values were used 



The Lambert scale A at nt = 5 is A = 0.2642 GeV in order for the expansion to reproduce the world average 
value a p t(M§ 
in our evaluations 

On the other hand, in the MA model [^-[3, the value Ama = 0.260 GeV (at n/ = 5) is the one that reproduces the 
high energy QCD phenomenology (see also [25|]). This value of Ama corresponds here to the Lambert scale value in 
MA Ama = (0. 260/0. 213)A = 0.3225 GeV. In MA, the renormalization scale invariant mass ml can be obtained by 
replacing in equation (p3"|) the noninteger powers a(m 2 y° +3 (where j — 0, ... ,4) by A^+jijn 2 ), the latter calculated 
via relations (|50p and (|10p using for the value of equation (|64p and for the discontinuity function the perturbative 
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expression p^' t \a) with the Lambert scale value Ama = 0.3225 GeV instead of A 
the case of MA in a value of rh^ somewhat lower than the one given in equation 

m[ MA) = 15.029 GeV . 



0.2642 GeV. This then results in 

(73) 



On the other hand, the usual perturbative QCD (pQCD) approach in evaluating the mentioned decay width is 
obtained by using the pQCD expansion ([58)) in powers of a(s) (p 2 = s = Mjj), with the overall factor fn 2 (s) there 
given by equation (|63|) . and the coefficients t n obtained from coefficients fi, "/j and Mj by using the integral relation 
([59")) . On both sides of equation (|59")) expansions in powers of a(p 2 ) at the fixed renormalization scale p 2 — Q 2 are 
used, and m 2 (s)/rfi 2 (Q 2 ) is also expanded in powers of a(Q 2 ). This then involves integrations of powers of (large) 
logarithms t = ln(s/Q 2 ) 



In = 



ds In™ (s/Q 2 



Q 



2)2 



(74) 



which are: I2 = tt 2 /3, I4 = 77r 4 /15, etc.; hk+i = 0). For details, see 
relations between i„'s and fk's and 7j's [their equations (22)-(24)]. At rif 
(f n ,t n ) = (5.66667,5.66667); (42.032,29.1467); (353.229,41.7576); (3512.2, 



38( and |26( (App. A there), and their 



= 5, the coefficients t n and /„ compare: 
825.747); for n = 1,2,3,4, respectively. 
Here we see that the effects of l^k integrals tend to decrease the absolute values of \t n \ in comparison to /„ for n < 4, 
and this makes the pQCD evaluation (|58[) numerically very well behaved. However, there appears to exist no reason 
for this tendency to persist at higher orders. Further, looking at the integrals (|74|) . we see that they involve integration 
over large RGE logarithms. 

The described pQCD method, i.e., the power series for T(a) of equation (|5"5|) truncated at j = 4, can be derived 



alternatively in the following way. We reorganize the power series for F(Q 2 ) of equation ([60]) into the form (|66[) in 
powers a pt {Q 2 )" 0+ i and include the terms up to j — 4. Then we put the resulting (truncated) series of F(Q 2 ) into 
the integral (|0 1 [) for T(er) where the integration path in the complex Q 2 plane is taken along the circular contour 
of radius a (the contour C2 in fig. 13}. However, instead of integrating the powers a pt (Q 2 ) Vo+: > as they are [with 



;-CT+l£ 



-a- ie 



Q -plane 



At: 



) 



FIG. 4: Paths Ci and C2 in the complex Q 2 -plane for the case when the coupling has unphysical (Landau) cut. 



Q 2 = a exp(i(j>) running], we use the perturbative RGE expansion of these powers around the fixed scale Q 2 = a > 0. 
This sometimes may not be the best approach, because the perturbative RGE expansion involves powers of relatively 
large logarithms ln(Q 2 jo) — i<fi, with — tt < <f> < ir. We thus obtain a series in powers a v t{<j) Va+: > which we truncate 
at j = 4. Reorganizing this series into the form with the overall factor m 2 (cr), we then obtain the series for T(a) of 
the form of equation (|58[) truncated at n = 4, where the coefficients t n turn out to be the aforementioned expressions 
involving /^'s and 7j's. 

The just mentioned pQCD method of equation (|58[) involves powers of a p t(p 2 ) at a fixed positive squared scale: 
p? = a = s = Mjj. On the other hand, our approach, equation (l3"5j) or (|54p [with T there defined via equation 

(|68p ], uses systematically the timelike quantities 2t^+i which are contour-integrated couplings A v +\ [cf. equation 



The value of rh b given in equation H65II is the one in pQCD and, to a large degree of precision, in all such analytic QCD models where 
the spacelike analytic coupling merges fast with a p t(Q 2 ) at high Q 2 : Ai(Q 2 ) — a p t(Q 2 ) ~ (A 2 /Q 2 ) n at Q 2 3> A 2 , with n > 3. One 
such model was constructed in [52ll . another in [53ll . both with n = 3. We note that in MA n = 1. 
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(|36|) ]. the latter corresponding to the generalization of the logarithmic "fractional" (noninteger) derivatives A n +i 
of the analytic coupling A\ [equation ([6])] in any analytic QCD model, or even in pQCD [equation (|25|) ]. We call 
our approch "fractional analytic approach" (FAA). It can be applied to evaluation of physical quantities [spacelike 
quantities T{Q 2 ), or timelike quantities T(s)] in any analytic QCD model. Further, it can be applied formally even to 
evaluation of high energy timelike quantities T(s) in pQCD, provided that q 2 — s is large enough: s > A 2 where A 2 
is the highest positive value of the nonphysical (Landau) cut of a pt (Q 2 ) in the complex Q 2 -plane. This is so because 
the contour integration in fig. [2] in sec. IIIII for equation (|36j) can also be applied to the pQCD coupling 

1 t—v+ie iq'2 

aw+iCa) = ^ J a i£ 2 ) , (75) 

provided the integration along the cut avoids the cut (including its unphysical part), see the modified path C\ in fig. [J 
(in comparison to the C\ in fig. [2]); and provided that at the same time the integration along the circular path C2 
gives the same result - the latter is the case only if a > A 2 . The numerical results, for T(s) and T(H — > bb)(s), as 
a function of the Higgs mass Mg = y/s, arc presented in figs. [5] a, b, respectively. We also include the curve for the 









0.045 


; s x (a): 7~ vs M„ 


- FAA L - 


0.044 




pQCD 1 ; 
MA 1 








0.043 


: FAA%. ' - - ^ 




0.042 


^^.pQCD "-MA 




0.041 






0.040 






0.039 







140 160 

ikf H (GeV) 



(b): r(MeV) 



pQcu 



FAA 

pQCD 

MA 



BAA 



M«(GeV) 



.''MA 



FIG. 5: (a) The (dimensionless) quantity 7~(s), as denned in equation J68I) . as a function of the Higgs mass Mh = for our approach 
(FAA) of equation (135 b . for the pQCD approach of equation 1581 (in both cases A = 0.213 GeV at rif = 5); and for the minimal analytic 
(MA) model with FAA approach (Ama = 0.260 GeV at n/ = 5); (b) the same but now for the decay width F(H — > bb). 



Minimal Analytic (MA) model, evaluated by our method and using A(n/ = 5) = 0.260 GeV. Our curve (FAA) can 
be interpreted as the result of application of our method in perturbative QCD, or in any such analytic QCD in which 
the values of the discontinuity function p\(cr) = A\{— a — it) do not differ from the values of the pQCD discontinuity 

function pf^ip) = a pt (—<J — ie) at yfa > Mh- 

We see that the FAA and pQCD curves are close to each other. The MA curve for T(H — > bb) comes close to 
the FAA and pQCD curves because the effects of different values of scales A (0.260 GeV instead of 0.213 GeV) and 
different values of rh 2 [equations ([73| instead of equation (|64|) ] tend to cancel each other. The evaluation for the 
curves was performed by using the renormalization scales such that |/i 2 | = s (= Mfj). If we vary the value of the 
renormalization scale |/i 2 | from s/2 to 2s around s, the values remain very stable; for example, when Mh = 150 GeV, 
the values of T are (4.062 ± 0.003) • 10~ 2 in the FAA method case, and (4.070 ± 0.001) • 10~ 2 in the pQCD method 
case. 

The decay width, by FAA method, is T(H -> bb) = 2.03, 2.82, 3.57 MeV, when M H = 100, 150, 200 GeV, respectively. 

It turns out that application of the expansion (|54|) in our FAA approach, instead of the expansion (|35"|) , gives the 
same result. Furthermore, the spacelike couplings A„(Q 2 ) obtained via equations ([501) and (fTU)) . and the timelike 
couplings 2l„ (Q 2 ) obtained via equations (|4"0"]) and (|4"4"]) and (fB"3"]) . when evaluated in such analytic QCD models which 
have pi(c) = Ima pt (— a — ie), i.e., in MA- type models, give results numerically indistinguishable from the expressions 
[see also equation 

Ai^HQ 2 ) EE - 

7T 

1 r°° r/rr 1 

ar T) M = - / ^a$ t (-cr'-ie), (77) 

7T J a <T 



da lmflp t (-(T 



(76) 



for any positive v (in general noninteger); and Q 2 in the complex plane outside the negative semiaxis: and a > 0. 
This is another check of consistency of our method, because it shows that the APT construction, 0, Q , of higher 
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power analogs in MA of Shirkov and Solovtsov, when generalized from integer to noninteger powers (n <-> v) as in 
equations (fT6" |) - (|T7)) . gives the same result as our method. 

Nonetheless, in conclusion, we stress that our method of construction of spacelike couplings A v and A v [equations 
(fTU)) and ([50])]. and of the corresponding timelike couplings 2t„ and % v [equations P0j) . (PPtj) . ([55)1 ] can be applied 
also to any other analytic QCD models, e.g., models where p\(p) ^ Ima p t(— cr — ie). The latter inequality can be 
expected in general at low positive a values, cf. [53 |. In such models, the APT- type construction, equations ([76"l) - ([77lh 
or modifications thereof, cannot be applied. 

We recall that the timelike couplings Silver) and % v (a) of our (FAA) method depend only on p\{cr') at cr' = cre w > a 

- see equations (HP)) , (UU, J53J). In the presented application of our method to T(H —> bb), however, a = s = > 
100 2 GeV 2 , i.e., only those Pi(cr') contribute for which cr' is very high (> Mfj). At such high cr' we can expect that 
pi(cr') = Ima pt (-(j' — ie). Therefore, in such cases the formulas (|77| for 21^'s can be applied [but not the formulas 
(TfBl for A^s] and they give the same result for 21^ (cr) as our approach. Really, the FAA curve and the MA curve; 
in figs.[5]are numerically reproduced by application of equations (|77|) . using the corresponding values A = 0.213 GeV 
and 0.260 GeV, respectively. 

It would also be interesting to apply our method to evaluation of low-energy timelike observables, where pi{cr) may 
differ significantly from the perturbative value. The method can also be applied to evaluation of spacelike quantities. 

V. CONCLUSIONS 

We presented a method of calculating spacelike and timelike QCD observables whose perturbation expansion in 
perturbative QCD (pQCD) has noninteger powers of the perturbative coupling a pt (= a s /n). 

The method can be applied in any analytic QCD model, i.e., (a) in any model with a given analytic spacelike 
coupling A\(Q 2 ) (where A\(Q 2 ) is the analytic analog of spacelike 20 a pt (Q 2 )); (b) or with a given discontinuity 
function pi(a) = ImAi(—cr — ie) (where a > 0). Specifically, first we constructed the analytic analogs A u +i(Q 2 ) 
of the (^-noninteger extension) of the logarithmic derivatives a p t,v+i{Q 2 ) oc d l/ a pt (Q 2 )/d(\nQ 2 )' J+1 , where v can be 
any real number larger than —2, cf. equations (|10l) . (|22|) . Furthermore, we constructed the corresponding timclikc 
(Minkowskian) couplings 2t I /+i(cr) (cr > 0), cf. equations (|43|). (J49j) . Subsequently, we obtained the analytic spacelike 
couplings Av(Q 2 ) as a linear combination of the aforementioned -4 y + m (Q 2 )'s (m = 0, 1, 2, . . .), where the couplings 
A U (Q 2 ) are analytic analogs (in any given analytic QCD models) of the powers a pt (Q 2 ) I/ (y any real number above 

— 1), cf. equation (f50|) . Furthermore, the corresponding timelike (Minkowskian) power analogs 2t I/ (cr) were constructed, 
as the corresponding linear combination of 2l 1/ + m (<5 2 ) (m = 0, 1, 2, . . .), cf. 

We further demonstrated that in the Minimal Analytic model (MA, also named APT) of Shirkov, Solovtsov, 
Solovtsova and Milton bMp , our method gives the same explicit results for A u +i(Q 2 ) and 2l„+i(er) at the one-loop 
level as the method of [24 - 1261 ] of Bakulev, Mikhailov and Stefanis (BMS; whose method can be applied in these 
MA- type models only), cf. equations (|17|) and (|47|) . When going beyond the one-loop level within MA, the explicit 
formulas for A y +\{Q 2 ) and 2l„+i (a) in j23 - [26| become complicated and the comparison with our results becomes 
harder. Numerically, though, we have strong indications that within MA both BMS and our method agree also 
beyond the one- loop level. We recall that our results are given in form of integral, i.e., they are less explicit than the 
results of [2^ - [26j for MA. However, our results are applicable in any analytic QCD, and look simple in its (integral) 
form. 

When the analytic QCD model is based on a beta function (3(Ai) which is analytic at Ai = ( 03 ) . i.e., in 
perturbative analytic QCD, we simply obtain A n = A™ [l|| and A v — A\. 

Furthermore, our method can be applied to evaluation of timelike observables T(s) within the nonanalytic pQCD, 
provided that s > A^, where A^ is the positive branching point of the unphysical (Landau) Q 2 -cut (0, A 2 ) of a p t(Q 2 ). 
This latter approach can be described as a RGE-resummed contour method, as opposed to the more usual fixed-scale 
contour method in pQCD. On the other hand, for spacelike observables, while our method can be applied in (any) 
analytic QCD, it cannot be applied in nonanalytic pQCD, and the results are different from the usual (nonanalytic) 
pQCD evaluation results. 

Further, if we work in an analytic QCD model for which the discontinuity function pi(cr) is equal to its pQCD 

counterpart p i ^ t \a) for large enough a > Mq (where Mq ~ 1 GeV > Al, Mo being a typical scale of the onset 
of pQCD), then the method gives for timelike observables T(s) at s > Mq the same result as the method gives in 



spacelike in the sense that Q 2 (= —q 2 ) G C\(— oo,0] 
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nonanalytic pQCD. For spacelike observables no analogous statement holds. This is so because the spacelike couplings 
[Ai(Q 2 ), A„(Q 2 ), A l/ {Q 2 )} are represented by integrals involving the values of pi{cr) along the entire cut of Ai(Q 2 ), 
while the timelike couplings [2li(s), 2l„(s), 2lj,(s)] involve only pi(a) for the cut sector a € (s,+oo). 

We applied the method to evaluation of the Higgs decay width into bb pair T(H —> bb), as a function of the Higgs mass 

Mh- The results of this evaluation turn out to be the same in pQCD and in any analytic QCD with p\{a) — Pi\cr) 
at a > Mjj, because this is a high-energy timelike observable: T(H — > bb) oc T(s) with s = Mfj ^> A 2 . 

It would be also interesting to apply our method in analytic QCD models to evaluation of low-energy timelike 
observables T(s) that involve nonintegcr powers, where pi(c) at a ~ s may differ significantly from the perturbative 
value; and to evaluation of (low-energy) spacelike quantities in such models. For the latter, we plan to investigate the 
structure functions of the deep inelastic lepton-hadron scattering in analytic QCD models. 
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Appendix A: Coefficients k m (v) and k m iy) 
1. Coefficients kmiy): results 

We need to find the coefficients k m (i>), where m = 1, 2, . . ., and v is any real number. These coefficients are derived 
later in this appendix. First we write down their solution explicitly. 

It turns out that they involve derivatives Z m (y + n) of the Riemannian T-functions 

ZM s WTT)d^\ r(i + x) ) ,=o ' (A1) 

These functions can be expressed in the form of the Euler ^-functions and their derivatives 

d d m 
ee ^-hvT(v) , ^ {m \v) = -7— *M . (A2) 
w dv w w dv m w y ' 

S x {v) = , (A3) 

(_-\\m—l , , 

S m iy) = \ 7Tj (^ (m - 1) (^ + l)-^ (m - 1) (l)) (m = 2,3,...). (A4) 
(m — 1)! v / 

We note that S m (y) for integer v — n (and m = 1,2, . . .) coincide with the usual harmonic numbers (of order m): 21 
S m (n) = J2k=i k~ m - I n terms of these functions, the functions Z m (v) of equation (|A1[) are 

Zo(iz) = 1 , Z x {v) = Si(v) , Z 2 (v) = Sxiyf - S 2 {v) , 
ZM = S 1 (v) 3 -3S 1 (v)S 2 (u) + 2S 3 (v) , 

Z A (y) = 5i(^) 4 - QS l {v) 2 S 2 {v) + iS 2 {vf + 85i(i/)S 3 (i/) - 6S A (u) , (A5) 
In terms of the quantities Z m (y), as given by equation (|A1[) . the coefficients k m [y) are 

ki(v) = ucxB^u), (A6) 

k 2 {v) = v{v + l)(c 2 B 2 (v) + ^B ltl (v)\, 
Hv) = V{V + ^ + 2) U BM + Cl c 2 B lt2 {y) + | B 1>ltl (u)j , (A8) 

W = v{v + 1)(l/ 6 + 2)(t/ + 3) U BM + 4 b.,M + °-f b x ,m + ^ Bi^v) + 1 iW")) , (A9) 



(A7) 



In Mathematica [28H . 9jf( m )(i/) is denoted PolyGamma[m, v]; and S m {v) is denoted HarmonicNumber[n, m]. 
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where 



Bi(v 
Bi,x,2{v 



Z x {v)-1, B 2 {y) = 



v - 1 
2(i/ + 1) 



B 1<1 (p) = Z 2 (u+l)-2Z 1 (u) + l , 



11 1 v , N 2 1 11 

6 i/ + 1 ^ + 2 i/ + 2 i/ + 1 v + 2 6 



= Z 3 (^ + 2)-3Z 2 (^ + l) + 3Z 1 (^)-l 



1 2 

12 ~ („ + l)(„ + 2)( I / + 3) 
6 



r, / N 5 1 



1 



= 1- 



(^ + 2)(^ + 3) 
3(i/ +1) 



^ + 2)+ — — + 



+ + 

12 



12 + l i/ + 2 i/ + 3 ' 
4 1 1 13 



z/+l ^ + 2 i/ + 3 6 

6 \ „ , . 6 



11 + 



11 38 



^ + 3 \^+2 z/ + 3 

= Z 4 (i/ + 3)-4Z 3 (i/ + 2) + 6Z 2 (i' + l) -4Z x (v) + l 



v+1 v + 2 v+Z 3' 



(A10) 



2. Coefficients k m (v): results 

Inversion of the relations (f28| gives us the expansion (|29|) and the k m (v) coefficients in terms of the previously 
written coefficients ki(y + n) (£ = 1, . . . , m, n = 0, 1, . . . , m — 1) 

fei(z/) = , fc 2 (^) = k x {v)k x {v + 1) - fca(i/) , 

= -k x {v)k x (v + l)h(u + 2) + k x (v)k 2 (v + l) + ki(u + 2)k 2 {v) - k 3 (v) , 
k±{v) = fciO)fci(v + + 2)fc 1 (i/ + 3) - k x {v)k x {v + l)k 2 {v + 2) - k x {v)k x {v + 3)k 2 {v + 1) 

-k x {v + 2)ki(v + 2>)k 2 {v) + k 2 (y)k 2 {v + 2) + k x (y)k 3 {v + 1) + k x {v + 3)fcs(i/) - fc 4 (f) • (All) 

It is possible to check that our formulas, equations (|A6|) - (|A9[) and (jAllj) . give 

fc n (l) = k n (l) = , (A12) 

reflecting the fact that, by definition [equations (f2"5]) and ©J, _4i = A x and a pt ,i = a pt . 

It is interesting that the coefficients k m (y) can be cast in an equivalent alternative form which is somewhat similar 
to the expressions for the coefficients k m {v), equations (|A6[) - (|A10[) . and where instead of the derivatives (|A1|) . another 
type of derivatives appears naturally: 



T(l-x) 



x=Q 



(A13) 



dx m \T{v + l-x) 

In terms of the harmonic numbers of order k, Sk(n) = Y^— x s~ k , functions Z m (n) and Z m (n) can be expressed as 



Z m {n) 



1 d m (T{n + l + x) 



Z m (n) = T(ri + 1) 



T(n + 1) dx m V T(l + x) 
d m ( T(l - x) 



dx m \T(n+ 1 - x) 



x=o dx r 
d m 



■ exp 



E 

,fe=i 



(-l) fe+1 x fe 



S k (n) 



x=0 



x=o dx r 



■ exp 



,fc=i 



(A14) 
(A15) 



This means that for Z m 's, formulas very similar to those of equations (|A5|) for Z m 's, are valid 

Z {v) = 1(= Z {v)) , = ft (i/)(= Z x {v)) , Z 2 (v) = S x {vf + S 2 (v) , 

ZM = S x (uf+3S 1 (u)S 2 (u)+2S 3 (u) , 

Z 4 (i/) = 5i(^ 4 + 65iH 2 5 2 (z/)+35 2 (^) 2 +85iH5 3 (^)+654(^) , 



(A16) 
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The coefficients k m {v) are then written in a form very similar to the formulas (|A6[) - (|A9[) for fc m (f)'s 

k x (v) = -vciB x {v), (A17) 
k 2 (v) = u{v + l)(^-c 2 B 2 (u) + ^-B 1>1 (u)^, (A18) 

Hv) = viy + l) ^ V + 2) ^ - c 3 B 3 {v) + cic 2 Bi,a(i/) - ^ . ( Al 9) 

W = Kiy + 1)( 7 2)(iy + 3) ( - c 4 i^) + c 2 S 2i2 (,) + ^ 5 llS M-^ S M , 2 (-) + | Bi,i,i,iM) , (A20) 

B x [y) = Z x (v)-\ , B 2 {v) = B 1 , 1 (y)=Z 2 (v)-2Z 1 (v + l) + l, 

2[y + 1) 

Mv) = \- — l-r + -^r , B 1<2 {v) = Uzxiv + 1) - 1 1 



where 



6 v + 1 v + 2' , w 6(z/ + l) V v ' v + 2 

= Z 3 {y)-ZZ 2 {y + \)+$Zi(u + 2)-\ , 

~ , s 1 2 ~ , , 13 1 1 1 

B&) = — - . ■ ■ ■ - , B 2 , 2 (u) = — - 



12 (v + l)(i/ + 2)(i/ + 3) ' w 12 i/+l ^ + 2 i/ + 3 

Bi.afr) = fl- 7 , , 5^ Z 1 (v + 3) + l+ 4 



(z/ + l)(i/ + 2)y v ' 6 i/+l ^ + 2 i/ + 3 ' 
Si^M = ^f^* + 2) - ( 1 - , , J Z x {y + + 2 ' 1 



v+l v 7 V + + 3 i/ + 2 i/ + 3 ' 

Bi.i.i.iM = ^M-4Z 3 ^ + l) + 62 2 ^ + 2)-4Z 1 (i/+3) + l. (A21) 

3. Reorganization of the power series; RS-dependence of new coefficients 

The relation (|29p between a^ t and a p t,!/+m's allows us to obtain immediately the relation between the coefficients 



Fj of the usual perturbative expansion (|23|) and the coefficients of the reorganized ( "modified" ) expansion (f30l) 

= 7i+fci(z/ ) , (A22) 

F 2 = F 2 + fci(z/ + 1) J"i + fc 2 M , (A23) 

F 3 = ^3 + *i(«o + 2)J" 2 + fc 2 (^o + + fc 3 (fo) , (A24) 

Fa = Fa + h(u + 3i)F 3 + k 2 (u + 2)F 2 + h(v + l)Fi + fc 4 K) , (A25) 

All these are coefficients when the renormalization scale (RScl) /i 2 is taken to be /i 2 = Q 2 . If using another RScl 
M 2 Q 2 to calculate the observable F(Q 2 ), we can obtain the RScl dependence of the coefficients Fj(/j, 2 /Q 2 ) by using 
the relations dA v {n 2 ) / 'din /i 2 = — (3qvA u +\ ([i 2 ) [or equivalently: (ia p t,i/(/i 2 )/dln/i 2 = — /3o^a p t,„+i(/i 2 )], cf. equation 
ITl]) . and the RScl independence of the (spacelike observable) F(Q 2 ). The resulting expressions are 

ft&W) - A + t^^^S I"' (|) *U ■ (A26) 

where we denote throughout Fk = 

Since the coefficients fci(^o + are RScl independent, the relations (|A22I) - (|A25|) are valid at any RScl /i 2 , i.e., we 
can replace there Fj n- Fj(fi 2 /Q 2 ) and Fj n- Fj(p? /Q 2 ). 

4. Derivation of expressions for k m (v) 



In this Subsection, we will omit the subscript pt (a p t <-> a). Let's consider the relation 

(-1) d 
n/?o d (In Q 



a n+1 (Q 2 ) = { — ^ ) H„(Q 2 ), (A27) 



which is a consequence of the definition (|25l) . In addition, we will use the following property: 

(a - Tn ^2 J a " 1 = ma m+1 (l + c x a + c 2 a 2 + c 3 a 3 + c 4 a 4 + . . . ) , 
V Pa d{lnQ 2 )J 

which follows from the RGE, equation (f2l)|) . We then obtain 

^ ,,* 2 . ) a n (Q 2 ) = a n+1 (l + Ki(n)o 1 +K 3 (n)o a +X:,(n)a 8 + C 1 (n)fl 4 

■ (l + cia + c 2 a 2 + c 3 a 3 + c 4 a 4 + . . . ) , 



where 



Equation (|A28[) can be rewritten as 

a„+i(Q 2 ) = a ,l+1 ^l + [/Ci(n)+ci]a 1 + [/C 2 (n)+ci/Ci(n)+c 2 ]a : 
+ [£3(71) + cxK.2(n) + c 2 K,\(n) + C3] a 3 
+ [/C 4 (n) + ci/C3(n) + c 2 K, 2 (n) + C3/Ci(n) + C4] a 3 + . . . 

Comparing equations (|28|) and (|A31|) we obtain the following recursive relations: 

(n + 1) 

fci(n + l) = -A;i(n) + ci, 

n 

(n + 2) , , , (n + I) , , „ 

fc 2 (n+l) = -fc 2 (n) + ci- -fci(n) + c 2 , 

rt n 

1 / ^ (« + 3), , , (ra + 2), . . (n + 1), , . 

k 3 (n + l) = - -fc 3 (n) + ci- -£ 2 (n) + c 2 ^ ^fci(n)+c 3 , 

n n n 

, , ^ (« + 4 ), / x fa + 3), , N (n + 2). , , (n + 1), , , 

fc 4 (rc+l) = ^ ^fc 4 (n) + ci^ ^/c 3 (n) + c 2 ^ -fc 2 fa)+C3" -fcifa) + c 4 . 

n n n n 

We note that 

kj(2) = Cj , fcj(l) = 0. 

In order to solve the recursion relations (|A32[) - (j A35|l . we find first solution to the following general recursion 

n + ai 

k{n + 1) = k(n) + c(n) , 

n + a.% 

where c(n) is some function of n and a% and a 2 are some parameters. 
If c(n) = 0, then the solution of equation (|A37[) is very simple 

1 (n + a 2 ) 

where c is an arbitrary constant and T is the Riemannian T-function. 

If c(n) 7^ 0, it is convenient to introduce a new variable k(n) which is related with k(n) as 

T(n + a 2 ) 

Using equation (IA39[) in equation (|A37I) we obtain 

L(n + ai + l) 



21 



with the solution 

if \ l( \ j V r(j + a 2 + 1) 

where s is a chosen number. Below it will be convenient to use s = 2, because fcj(2) = Cj. 
So, for r(n), we have 



(A41) 



fc(n) = 



T(n + ai) 
T(n + a 2 ) 

T(n + «i) 
r(n + a 2 ) 



r( s + a 2 ) 

r(s + ai) 

T(2 + a 2 ) 
T(2 + ai) 



k( S )+ J2 



j=s+ 



j=3 



^4^-1) 

i r(j +ai) 
Sf^4c (j - 1) 



(A42) 



Having the solution (|A42[) to the recursion relation (|A37[) . we can proceed to solving the recursion relations (|A32[1 - 
(lAl5l) . 

1. Consider the recursion (|A32|) : it corresponds to the general case with 



k = ki, a\ = 1, a 2 = 0, c(n) ~ ci, fci(2) = c\ . 
So, using the solution (|A42[) . and fci(2) = ci, we have 



k\(n) = nci 



where 



n 1 

.(«) = F- 



(A43) 



(A44) 



(A45) 



are harmonic numbers (of m'th order), which can be related with the (to — l)'th derivative of ^-function: 

I _ 1 \m , v jm 

S m+1 (n) = ±-L- (^ m \n + l)-^ m \l)), vpM( x + i) = -^(z + 1) . 
m! v / ax 

The ^-function is in turn the logarithmic derivative of the corresponding T-function: 

*(x + l) = — lnlYx + l) 

ax 

and 7 = — is Euler constant. 

So, equation (IA44[) can be represented in the following form: 



(A46) 



(A47) 



ki{n) = nci *(n + l)-tf(l)- 1 = nci *(n + 1) - *(2) , 



(A48) 



which is well-defined also for noninteger values n i— > v 



k x {v) = vc x ¥(i/ + l)-¥(2) , 



(A49) 



It can be considered as the analytic continuation of the coefficients k m (n) (-> k m (u) based on the corresponding procedure for harmonic 
numbers |57|1 . 
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2. Consider the recursion (|A33|) : it corresponds to the general case with 

eti = 2, a 2 = 0, fc 2 (2)=c 2 , 
n + 1 

c(n) = ci • fci(n) + c 2 = c? • (n + l)[5i(n) - 1] + c 2 

n 



So, using the solution (IA42I) . we have 
The coefficient in front of c 2 has the form 



(») = n(n+l) ^2+g -T^ T y |c 2 +^ 



jfrtf-l)-!] 



6 + §G J' + l) 6 + ( 



n n+1 



E 53 ) j 6 + (3 n + 1 



3=3 J=4 . 



n- 1 
2(n + l) 



The coefhcient in front of c 2 has the form 



Here, the second term on the right-hand side is 

n 

Y,~ = Si(n)-5 X (2) = S^n) 



i=3 



3 

2 ' 



while the first term is 



^ 7 + 1 ^ 



J'=3 



j=4 



^(n + l) - S 2 (n + 1) J - 1 = -Z 2 (n + 1) - 1, 



where 



Z m (n) 



1 d' m /T(ra+l+x) 



a=o ciir r 



■ exp 



.k=l 



T(n+1) (dx m \ T(l + x) 

For several first values, Z m (n) are 

Z (n) = 1, Zi(n) = Si(n), Z 2 (n) = Si(n) 2 - S 2 (n) , 
Z a (n) = S^n) 3 - SS^S^n) + 2S 3 (n) , 

Z 4 (n) = 5i(n) 4 - 6S*i(n) 2 S* 2 (n) + 8Si(n)S a (n) + 3S' 2 (n) 2 - 65 4 (n) 

In the case noninteger values n 1— > f, the (m'th) order harmonic numbers (|A45[) arc 

(_l)m-i 



■5 fc (n) 



S x (i/) = + S TO (i/) = ) j + 

(m — 1)! V 



where 



*(!/) 



— in(r M) , = _*(„) 



So, the result for fc 2 (n) has the form 



fc 2 (n) 



n(n — 1) n(n + 1) 



c 2 



2 2 
which is well-defined also for noninteger values n n- v 



v{v-l) v(v + l) 



■ c 2 



c?-(Z 2 (n + l)-2Zi(n) + l) 



^•(Z 2 (f + l)-2Zi(i/) + l) 



(A50) 



(A51) 



(A52) 



(A53) 



(A54) 



(A55) 



(A56) 



(A57) 
(A58) 
(A59) 

(A60) 
(A61) 



3. The results for ^3(71) and ki(n) can be obtained similarly. After the replacement n <— > v they have the forms 
as given in equations ([A"8" ]) -([A"9 )) and (|A10j) . 
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Appendix B: Large- behavior of coefficients 



Here we calculate the asymptotical results for the coefficients T n at n — > oo of the expansion (|30p . assuming the 
standard Lipatov-type behavior 34] J- n ~ n! for the coefficients T n of the original expansion 
It is convenient to use the following form for the coefficients T n at n — > oo: 



_ r(n + z/ ) 



where 6 ~ 1 and T is the Euler Gamma function. 
From equations (|A22I) - (|A25[) we conclude that 



T n = k n - m (u + m)F m , k = 1, J"o = 1 



m=0 



Firstly consider the first several terms on the right-hand side. 

Let in — Ti — 1. Due to equations (|A17I) and (IA2ip . the right-hand side of (|B2[) has the following form: 



ki(v +n— l)-F n -i = —c\v Z\{y) - 1 



i> ) 



-Ci&' 



It can be rewritten also as 

, Tn-\ 



k\{v§ + n - 1)'- 



£l 
' 6 



Zi(m + n-l)-l 



c\ ( d 



b \dx 



i> ) 
r(i-x)r(^ + n) 

r(^ + n — a;) 



Zi(i/ + n — 1) - 1 



a;=0 



where we use equation (|A13[) on the right-hand side. 
At m = n — 2 we have from ||aT8|) and (IMTj) 



kfa + n - 2)%^ = if-ca ^(^o + n - 2) + ^ Bi.ifa + n- 2) 



where at n — > oo (see equation (|A21| ) 



(Bl) 



(B2) 



B 2 (v + n - 2) w - , Bi,i(i/ + n- 2) « Z 2 (2/ + n - 1) - 2Zi(^ + n - 1) + 1 . 

We use the simbol w to show the asymptotics at n — > oo. In particular, in the B^i-case, the symbol rs involves the 
replacement of the argument + it - 2 in Z2 by i/o + n - l. 23 Similar replacements will be used below. 
So, we have 



^2(^0 + n — 2) 



1 



-F„ 26 2 

When m = n — 3, using (|A19[) and (|A21[) we have 



T{l-x)T(v + n) 
r(v + n — x) 



x=0 



k 3 (v +n-3)'- 



1 / ~ c? ~ 

^ ( - c 3 B 3 (^ + n - 3) + cic 2 B lfi (v + n - 3) - y £1,1, 1(2^0 + n - 3) 



where at 



(see equation (|A21|) ) 



#1,2(^0 + n- 3) w ( Zi(i/ +n- 1) - i ) 



-63(^0 +n- 3) 

^1,1,1(^0 + n — 3) « ^3(^0 + n - 1) - 3 Z 2 {va + n - 1) + 3 Zi(v a + n - 1) - 1 



It can be checked that Z^{n + m) / ' Z^(n) — > 1 when n — >• 00. 
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Hence we have 



fo(vo + n — 3) 



1 



C3 ^^ ( d r , 

y + 3ciC2 U-6 1 " 



Finally, when m = n — 4, equations (|A19I) and (1A21[) lead to 



dx 



- 1 



r(l - x)T{v + n) 
r(^ + n — x) 



x=0 



~ J- A 1 / ~ ~ 

fc 4 (^o + n - 4)-^- = g^i ( - c 4 #4^0 + n - 4) + c\ B 2 ,i{vo + n - 4) + £1,3(^0 + n - 4) 

2 4 

-YBi.i.j^o + n - 4) + ^ + n - 4) 



where at n — » oo (see equation (|A21[I ) 

1 



~ 13 ~ ~ 1 

^2 ' ^2,2(^0 +n-4) w — , Si )3 (^)+n-4)RiZi(vb + n-l) + - 



£4(^0 + «- 4) 

•§1,1,2(^0 + n - 4) « 3Z 2 (^ + n - 1) - Zi(f + n - 1) + | , 
Bi,i,i,i(v + n -4) w Z 4 (^ + n- 1) - 42 3 (i/ + n- 1) + 6 Z 2 (>o + n - 1) -4Zi(i/ + n- 1) + 1 
So, we have 

J r „_4 1 / 13d 2 — C4 



£4(^0 + n - 4) : 



246 4 



+c ^Tx~ x 



— - 

6 



2c i c 3 ( — + 7 ) -6c?c 2 



_d_ _ l 

da; 6 



31 
3G 



T{1 ~ x)T(v + n) 
T(uq + n — x) 



x=0 



(B3) 



(B4) 



1. Contributions of the powers of ci 

Taking only the terms ~ c™ (m = 1, n), we have 

T(l -x)T(yo +n) 



- T V 



(-i) m cr M _i 

m! 6 m \ c?x 



r(i/ + n — x) 



(B5) 



At the beginning it is convenient to consider the sum on the right-hand side at n — > cxd. Moreover, the last term can 
be represented as 

r(l-x)T(u 0+ n) = _ fdyy-° { l-yy^-\ 

r(z/ + n — x) J 

Application of the operator (d/dx — (3) to the term y~ x has the simple form 

dx J \ y 



a/b e ci/3/b 



Then we have for the series on the right-hand side of (|B5|) at n —> 00 



» (-irc f / i_ ■ 

^ m!6 m V y 

m=0 x H 

7(1) 



exp 



Thus, the contribution J-"„ has the following form after integration on y 

F {1) - T R 



l/6 r(i + ci/6)r(^ + n) 

T(i/ + n + ci/b) 



where the operation R n [F(ci)} takes the first (n + 1) terms of the expansion of F in powers of c\, and R n (x) —> x 
when n — > 00. 
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Since n is large, the difference between R n [F(ci)) and F is very small (~ l/(n + 1)0- So, we can omit the R n 
operation and take the contribution of the terms ~ c™ (m = 1, n) in the form 

^^^ r ^ +c ^°+;) . ( B6) 

r(z/ + n + ci/o) 

Note that the contribution JF^ is very important because the coefficient c\ is universal and has nonzero value in 
MS-like schemes. In the scheme where all Cj = {j > 2), Tn^ represents the full contribution to the reexpression of 
the expansion (|2"5|) to the one of equation (|30p : = T n , As we can see in the next subsection, the contributions 

~ c j (j > 2) can be expressed also in the form ~ . 
Use of Stirling's formula in (|B6|) gives us 

^/T n ^e^ b T{l + Cl /b)^. (B7) 

When the coefficients in the original series are nonalternating in sign, such as the one encountered in the Higgs decay 
width, equation (1671) . b is positive and the above ratio even tends to zero when n increases (we note that c\ is positive 
for all n/ < 6). 

Thus, the Lipatov-type asymptotics takes place for both T n and T n coefficients in the Cj = (J > 2) scheme: only 
the subasymptotical terms are changed. 



2. Contributions of ~ Cj (j > 2) 



The terms proportional to the first power of ci have the following form 

2 



C2 , CiC 2 



(- 




9 

CJC2 


(- 




\dx 


-I) 




\dx 


4) 
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36 



T(l-x)T(v + n) 
r(j/ + n — x) 



It is convenient at the first stage to exclude the term ~ (31/36)cfc2 from the consideration. It will be considered 
below together with the term ~ C4. 



Without the term ~ (31/36)c^C2, the contribution has the following form 



r(2) 



_C2_ 

2b 2 



(_l)m c m 
ml b m 



d__l 

dx 6 



T{l-x)T{u + n) 
T(v + n — x) 



x=0 



(B8) 



Repeating the calculations for J 7 ^ done in the previous subsection we obtain for 



2& 2 n n ~ 2 



ei/(6t) r(i + Cl /ft)r(m + n ) 
r(i/o + 71 + c\/b) 



where, as in the previous subsection, the operation R n _2[F(ci)] takes the first (n — 1) terms of the expansion of F on 
c x . 

Since n is large, the difference between R n -2[F (ci)] and F is very small. So, we can omit the R n -2 operation and 
take the contribution of the terms ~ c™ (m = 1, n — 2) in the form 



fi-2) „ Cl/ ( 66) r(l + Ci/6)r(l/ + n) ^ __^_TW p -5 Cl /(6b) 

■ r " ~ 2fo2 - r « e r ( jyo + n + Cl / 6 ) ~ 2b 2 " 



(B9) 



Taking the terms ~ C3, we have 



C1C3 



126 3 126 4 Vdx 6 



1 



C3 

126 3 



•Fn 1 



(-cx) 1 fd 1 



1! 6 1 V da; 6 



r(l - iQl> + n) 
r(^ + 77 - x) 
r(i - x)iVo + n) 
r(fo + 77 - x) 



z=0 



a:=0 



Repeating the above calcuations, we obtain 

rO) w c Jl_ t r 

126 3 



; _ ci/(6fe) r(i + Cl /b)i> + n) 

T(t'o + 77 + ci/6) 
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and, because n is large, 

r^n-^Pe- 7 ^™. (BIO) 

Now we consider the remaining terms ~ l/b 4 . At the leading order (in ci), they contribute in two places: as the 
term ~ (13c 2 , — C4) in equation (|B4|) and as the term ~ (31/36)c 2 c 2 in the first equation of this subsection. Taking 
them together we have, at the leading order in c\ 

=(4) (2c 4 -26c 2 + 31c 2 c 2 ) 
" " 144& 4 • A " ' 

Adding to this the terms of relative higher order in ci, in analogy with above calculations for J 7 ^ (j = 2, 3), gives 

£(4) (2c 4 - 26c| + 31c 2 c 2 ) {je+1)ci/b 

To find the exact value of the factor k we should calculate the term k^^o + n — 5)F n -5/ F n in analogy with (|B4|) . It 
needs in turn the calculaton of the coefficients k§ and, thus, one step more in the analysis in appendix Hovewer, 
looking carefully at the above calculations, we note that the coefficients in the exponents, in front of —c\/b, rise with 
the index j of J 7 ^ ■ For the corresponding coefficient is equal to 7/6 and we suggest that k in (|B11|) should be 
bigger. 

So, we have for the coefficient T n the following approximation at large n values: 

J RijM/l C2 e -5 C1 /(66) C 3 -7 C1 /(6b) (2c 4 - 26c 2 , + 31c 2 C 2 ) ( fc +1)ci/b \ , , 

We see that the corrections from ~ Cj (j > 2) have the same sign and are decreasing in magnitude. Indeed, for 6 = 1, 
and when rif — 5, we have 

6=1, ci = 1.2609, c 2 = 1.4748, c 3 = 9.8357, c 4 « 86. , 
and these corrections apparently have decreasing magnitudes 

?n ~ Fn ] (l - 0.258 - 0.188 - 1.3 x (0.283) ( ^ +1) - , (B13) 

where exp[-(fc + l)c x ] « (0.283)^ +1 ' < 0.23 if k > 1/6. 

Thus, the contributions of ~ Cj (j > 2) are rather small and can be expressed through the contribution of ~ c±. 



Appendix C: Proof of equation (|43[) 



We prove the formula of equation (|43[) by mathematical induction with respect to n = 0, 1, For n = 0, this is 

the formula of equation (|40l) which was proven in the text. Now suppose that the formula equation (|43l) is valid for a 
given n. We will show that then it must be valid also for n + 1. 

If the formula of equation (|43l) is valid for a given n, we can use it and the recursion relation (|42|) to obtain 24 



2la+n+ 2 (c) 



(-l)sin(7r((5 + n)) 
^ 2 05 + n+l)(<5 + n)/3 o 5+n+1 



dw 

S+n 



W 



sin(7r(<5 - 



i9pi((re u ') dpi(a) wd 2 pi(a) 

diner ~ 1! d(lner) 2 ~ 
' dwpi{ae w ) 1 



w"- 1 d"pi(cr) 



dw 
1)) 



7T 2 (5- 



1 



ur 



S+n+1 (fi + n y+'. 



(n- 1)! d(lncr)" 
r pi(cre £ ) 



1 



(CI) 
dpi (a) 



{5 + n)(S + n 



d 2 Pl (a) 



- l)e<5+»-i diner 
d k pi(a) 



(6 + n)(S + n-2)l\ e 5 + n ~ 2 d(\n<r) 2 
1 d n Pl (a) 



(6 + n)(S + n- k)(k - 1)! e 5 + n ~ k d{\na) k 



(5 + n)(S)(n- 1)! e s d(lna)' 



-} 



(C2) 



A somewhat similar procedure, in the context of one-loop fractional perturbation theory, was performed in [3l| (see appendix A there). 
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Here it is understood that e — > +0; and in the step from equation (jCll) to equation (|C2[) we performed integration by 
parts in the first term. 

Now we use in pi(ae e ) Taylor expansion in logarithm of the argument (lner + e) 



Pi 



n ± Jk 



e k d k pi{a) 
^ fc! (din a) k 



0(e 



n+l\ 



in the above expression (|C2|) . and obtain 

sin(7r((5 + n + 1)) 



2l(5+ra+2 (c) 



n 2 (5 + n + l)p° 



S+n+l 



dwp 1 {ae w ) pi(a) 



ur 



8+n+l 



n 1 

' ,(5+n-fc/ 



fc=l 



(5 + n) 



f 



f 



(6 + n)e s+n 

d k pi(a) 



k\ [6 + n-k){k-l)\ 



{dlnaf 



Using the identity 



f 



f 



kl (5 + n - k){k - 1)\ 



1 {S + n) 
k\(S + n-k) ' 



we can rewrite the expression (|C4|) as 

sin(7r(5 + n + 1)) 



(Die 1 



2l<5+n+2(c) 



7T 2 ((5 + n+l)^ 



5+n+l 



fc=l 



jfe! (dlno-) fc 



+ C(e 1 -' 5 ). (C6) 



(C3) 



(C4) 



(C5) 



Since < <5 < 1, we are now allowed to take the limit e — > +0 in the above integral (e 1 5 — > 0) and we conclude that 
the identity is valid also for n + 1. This concludes the proof of identity via mathematical induction. 



Appendix D: Coefficients of expansion of MS squared mass 

Integration of the RGE's (|26p and (p32"j) gives for the MS squared running mass the solution in the form of expansion 
63|) . with the coefficients Aij there being (note that 70 = 1) 

Mi = -J-( Cl - 7l ) , (Dl) 



M 2 = \m\-±- ((c 2 - 72) - Cl ( Cl - 71)) , (D2) 
1 2 

X 3 = -r^i+^i^2-rfl-((C3-Ta)-ci(pa-72) + (q-Pa)(ci-7i)) . ( D3 ) 



Mi = ^Mi-M 2 1 M 2 + ^Ml + M 1 M 3 

(( C4 _ T*) - ci(c 3 - 73) + (cf - c 2 )(c 2 - 72) + (-c? + 2cic 2 - c 3 )(ci - 71)) • (D4) 
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